CdV 2 04: A rare example of a collinear multiferroic spinel 
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By studying the dielectric properties ol the geometrically frustrated spinel CdV^CU, we observe 
ferroelectricity developing at the transition into the collinear antiferromagnetic ground state. In this 
multilerroic spinel, lerroelectricity is driven by local magnetostriction and not by the more common 
scenario of spiral magnetism. The experimental findings are corroborated by ab-initio calculations 
of the electric polarization and the underlying spin and orbital order. The results point towards a 
charge rearrangement due to dimerization, where electronic correlations and the proximity to the 
insulator-metal transition play an important role. 



PACS numbers: 71.45.Gm, 71.10.Ca, 71.10.- 
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Spinels form a big class of materials - probably as big as 
perovskites. There are many magnetic materials among 
them, with rich magnetic properties. In contrast to per- 
ovskites, however, there are practically no ferroelectrics 
in this class. Why spinels are such "bad actors" as to 
ferroelectricity, is actually not clear; the frustrated na- 
ture of their 23-site sublattice can possibly play some role 
[![. Only recently, the magnetically-driven ferroelectric- 
ity was found in some spinels with spiral magnetic struc- 
tures, the best known example being CoCr 2 04 Q, and 
there are reports of ferroelectric-like properties in ferro- 
magnetic semiconductors HgCr 2 S4 or CdCr 2 S4 Q and in 
charge-ordered magnetite [^-la] . 

Here, we report the discovery of magnetically-driven 
ferroelectricity in a ternary spinel with a collinear mag- 
netic structure CdV 2 04 - apparently the first such exam- 
ple. We found this multiferroic behavior experimentally, 
confirmed the presence of polarization and estimated its 
magnitude theoretically, using ab-initio calculations. 

Vanadium spinels AV 2 04 (A=Cd, Zn, Mg) recently at- 
tracted a lot of attention due to their magnetic structure 
and due to a possibility of orbital ordering. When de- 
creasing the temperature, all three materials experience 
a structural transition from cubic to tetragonal symmetry 
0, and at lower temperatures they also show magnetic 
ordering. The magnetic structure is antiferromagnetic in 
ab (or xy) chains, but is tt44 m xz ~ an d yz-chains (see 
Fig. la. 

Properties of related systems, with magnetic A-sites, 
e.g. Mn or Fe, are similar, but the magnetic structures 
are more complicated. Several suggestions were proposed 
to explain the structural and magnetic transitions of V- 
spinels by different types of orbital ordering. V 3+ ions 
have two electrons in triply-degenerate ti g orbitals, which 
in the tetragonal phase are split into the lower occupied 



rry-orbital, and doubly-degenerate (zx, yz)-orbitals for 
the remaining second electron. At least three pictures 
of orbital ordering of this remaining electron were sug- 
gested: alternation of xz- and xy-orbitals in z-direction 
Q; occupation of complex orbitals xz ± iyz Q, and 
tetramerization (occupation xz-xz-yz-yz in xz- and yz- 
chains) ■ The description of orbital order should take 
into account the interplay between electron correlation, 
spin-orbit coupling, and cooperative JT distortions [llj ]. 
Recently, a novel picture was suggested on the basis of 
ab-initio calculations 12|: almost no orbital ordering of 



this second electron, but strong alternation of V-V dis- 
tances in xz- and yz-direction, with ferromagnetic bonds 
becoming shorter. 

The magnetic structure, with Tt44- spin ordering along 
a chain, is reminiscent of the situation in E-type mangan- 
ites like HoMn0 3 [HQ and in CasCoMnOfipj]. which 
are both multiferroic. This would suggest that multifer- 
roicity can be found in V-spinels as well. This is indeed 
confirmed by our experimental and theoretical study. 

The single-phase, poly-crystalline samples of CdV 2 04 
were prepared by solid-state reaction in evacuated quartz 
ampoules as described in Simultaneously, also the 

samples of ZnV 2 C>4 and MgV 2 C>4 were made. Structural 
and magnetic measurements confirmed the known behav- 
ior: in our samples of CdV 2 C>4 a cubic-tetragonal tran- 
sition occurs at Tg w 95 K, and at Tjy ~ 33 K there 
appears antiferromagnetic ordering. Both transitions 
cause pronounced anomalies in the thermal expansion 
coefficient a(T), measured in a home-built capacitance 
dilatometer (see Fig. [T^). The measurements of elec- 
tric polarization were made in a conventional 4 He-flow 
magneto-cryostat (Oxford) by evaluating the integrated 
pyro-current recorded with a high-precision electrome- 
ter. Also the complex, frequency dependent dielectric 
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FIG. 1: (Color online) Thermal expansion alpha in the whole 
temperature range (a) and in the vicinity of the magnetic 
transition (al). The complex dielectric permittivity e* of 
CdV2C>4 is displayed as real part e' (b) and as dielectric loss 
e" (c) for different frequencies between 82 Hz and 0.3 MHz 
measured with a stimulus of ~1 V/mm. The data displayed 
in frame bl) was measured in the vicinity of the orbital order 
transition at 1 Hz in external fields of zero and 10 T. Frame 
b2) shows the small dielectric anomaly at the structural tran- 
sition below 100 K. 



response was measured employing a frequency-response 
analyzer (Novocontrol) with pseudo four-point probing 
for the linear measurements under small stimulus and 
with two-point probing for the non-linear measurements 
in electric driving- helds up to 220 V rms . For both pur- 
poses, silver paint contacts were applied to the platelike 
poly-crystalline pellets in sandwich geometry with a typ- 
ical electrode area of A w 10 mm 2 and a thickness of 
d w 1 mm. The uncertainty in the determination of the 
exact geometry together with additional (but constant) 
contributions of stray capacitances in two-point probing 
results in an uncertainty in the absolute values for electric 
polarization and permittivity of up to 20 %. 

Results of the dielectric measurements are shown in 
Figs. Q] and El The dielectric permittivity displays a 
small anomaly at the structural transition temperature 
Tg ~ 95 K (Fig. [TJd2) and a sharp jump-like anomaly at 
the magnetic transition Tjv ~ 33 K. The latter slightly 
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FIG. 2: (Color online) Electric polarization (a) measured via 
the pyro-current under cooling with a rate of 2 K/min in an 
electric poling field of 220 V/mm and under heating with the 
same rate in zero external field. Frame b) displays the lin- 
ear component of the dielectric permittivity measured with 
a stimulus of 220 V ms /mm. The corresponding (third or- 
der) nonlinear component \e n i\ normalized on the linear term 
is displayed in c); the "floor" of roughly 0.8 x 10 -4 can be 
attributed to the harmonic distortion of the high- voltage am- 
plifier. 



depends on an external magnetic field as demonstrated 
by the data for zero field and 10 T displayed in Fig. [Tpl: 
the anomaly shifts slightly to lower temperatures under 
applying a magnetic field. At the same time, the per- 
mittivity possesses a distinct frequency dependent con- 
tribution superimposed on the described features: a step 
in the e'(T) is accompanied by a peak in e"(T) shift- 
ing to lower temperatures with decreasing frequency. A 
very similar relaxational contribution along the c-axis 
was found in multiferroic rare earth manganites. There, 
it was ascribed either to the freezing of an overdamped 
polar lattice-mode or, alternatively, one could think of 
the relaxation of localized polarons at defect states 



17]. 



In any case these data demonstrate that the dielectric 
properties are to some extent influenced by the experi- 
mental time window. 

That the sharp dielectric anomaly at Tn is connected 
with the onset of spontaneous electric polarization is 
demonstrated in Fig. ®l displaying the electric polar- 
ization P(T). The difference between the cooling data 
detected in the presence of an electric poling field, com- 
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The explanation of the appearance of polarization in 
CCIV2O4 has some analogy with the mechanism proposed 
for E-type manganites ll| lij]. The distortion, in par- 
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FIG. 3: (Color online) Schematic view of a) the V-O ionic ar- 
rangements with spin magnetic structure (see black arrows) 
and short (long) V-V bond lengths highlighted by dashed 
(bold solid) lines; a') oxygen-induced dipoles (see ellipses) and 
long bond-length V-0 (see black lines) - symmetry inequiv- 
alent oxygens are denoted with different symbols/colors; b) 
calculated polarization according to DFT+t/ (left) and HSE 
schemes (right); V and O arrangement for c) tttt an d c') 
tt-14- s Pi n ordering along [101] or [Oil]; d) values of ai, @2, 
«2, Pi angles as a function of U e ff- 



pared to the zero-field heating data, may be explained 
by field induced contributions like reversible domain ori- 
entation or the above mentioned relaxational contribu- 
tion. The latter can be noticed as a broad loss peak 
in e"(T) above 10 K in the plot of the linear compo- 
nent of the complex permittivity, which were measured 
with a large stimulating electric field and is displayed in 
Fig. ^jp. The behavior of the nonlinear component of 
the permittivity e n i evaluated from the third harmonic 
of the dielectric response can be taken as an additional 
evidence for the onset of ferroelectricity (Fig. Wp). The 
onset of the third-order nonlinearity is followed by a de- 
cay towards lower temperatures which is related to the 
increase of the coercitivity: for temperatures below 25 K, 
the driving field of 220 V, ms /mm cannot switch ferro- 
electric domains anymore. It has to be noted that the 
second harmonic contribution (not shown) stays within 
the harmonic distortion of the setup at T/v, as expected 
for a symmetric polarization response M(H). 

All these data unequivocally demonstrate that there 
appears spontaneous polarization in CdV204 below T/v. 
The magnitude of polarization obtained in our samples 
(Fig. rjk) is P ~ 5 fxC/m 2 . One should expect that 
in a single-domain single crystal the polarization should 
be much larger. Similar measurements on poly-crystals 
of ZnV2C>4 and MgV2C>4 did not show ferroelectric be- 
havior: polarization is absent (within the accuracy of 
P < 0.1 fxC/m 2 ). We discuss possible reasons for such a 
different behavior below. 



ticular the shifts of oxygens, are different in the vicinity 
of ferromagnetic bonds and antifcrromagnetic ones. The 
main feature of the corresponding lattice distortions is 
the dimerization of V-V bonds 12| , but the shifts of the 
oxygens are also present, which ultimately lead to polar- 
ization. A specific feature of spinels, i.e. the existence 
of two oxygens "attached" to each V-V bond, compared 
to one oxygen in E-type manganites, makes the picture 
more complicated, but the main physical mechanism is 
very similar. 

In order to shed light into the origin of the ferroelec- 
tricity, we performed ab-initio calculations using the pro- 
jector augmented-wave (PAW) method [l9| with plane- 
wave basis sets as implemented in the VASP program 
[20l | . To describe the correlation effects of V-3cZ electrons, 
we used the DFT+U method 2l|, for different values of 
U e ff — U — Jh ranging between and 8 eV, and the 
Heyd-Scuseria-Ernzerhof (HSE) hybrid functional [22j ]. 
The HSE has been shown to provide an accurate treat- 
ment of solid state systems when the delicate balance 
between itinerancy and localization of correlated states 

ts an important role [23j , as in the system under study 

Starting from the tetragonal experimental crystal 
structure with I&x/amd centric symmetry [24{, we per- 
formed ionic relaxations for different values of U e f / , keep- 
ing the V magnetic moments ordered as ttJ-4- along the 
[101] and [011] directions (see Fig.|3^,). The relaxed struc- 
ture shows the formation of short (S) and long (L) bonds 
between ft an d ti moments, respectively, see also [12. 
The situation seems to be different in MgV2C>4 [25[. Our 
relaxed structure has a I4i polar space group. For large 
values of U e t / , the general trend towards formation of V- 
V dimers is reduced. Note that, imposing the ferromag- 
netic spin ordering, the structure retains the inversion 
symmetry, and no dimerization is found for any values of 
U e ff. Further inspection of the relaxed polar ionic struc- 
tures shows that the spin, orbital and lattice degrees of 
freedom are coupled: in plane xz,yz long bonds are ori- 
ented as reminiscent of a Kugel-Khomskii model [2H and 



of the experimental structure proposed for ZnV2 04 [27 1 
(see Fig.rSh). For the relaxed structure, we evaluated the 
ferroelectric polarization Pr p (see Fig. [3]d) according to 
the Berry phase theory [28|. We obtained a finite Pbp, 
directed along the c-axis, which is strongly related to the 
formation of a dimerized structure at intermediate val- 
ues of U e ff. The electronic charge redistribution due to 
dimer formation gives rise to differences between the BP 
and the point charge model (PCM) values of P. The red 
curve in Fig. [3}d) is the total polarization, ionic plus elec- 
tronic (BP) contribution, i.e. Pbp, the grey one is the 
estimate of the point charge model (using the ionic nom- 
inal valencies). The fact that they are different suggests 
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that covalency effects are relevant. 

We now discuss the mechanism behind ferroelectric- 
ity. First of all, let us recall that if we consider a FM 
(tttt) state with the suggested experimental IA\/amd 
symmetry [24l]. all V-O-V angles are equivalent along the 
chain (see Fig. [3b). Note, however, that V-0 distances 
are slightly different, due to the peculiar coordination 
of the spinel structure: each O is an "apical" one with 
respect to one V ion and a "planar" one with respect 
to the other neighboring V ion. This does not exclude 
a priori some degree of orbital-ordering, even in a FM 
spin configuration. As expected from the centrosymmet- 
ric space group, no polarization is found from our calcu- 
lations for this case. However, once the magnetic order 
ttii is imposed, (i) the angles ot\ and /3 2 («2 and f3\) be- 
come inequivalent as a consequence of the formation of 
short/long V-V bonds; (ii) «i and (3\ become different. 
The long V-0 bonds are arranged as shown in Fig. [3^'. 
This pattern is compatible with a weakly staggered xz,yz 
orbital ordering. As a result, the dipole moments P and 
P , originating from the inequi valency of oxygens, and 
schematically shown in Fig. appear due to different 
ai and («2 and (3\) angles; since P and P don't 
compensate, they give rise to a net P. This picture is 
valid for intermediate values of U e ff. At large values 
of U e //, however, the inequivalence of ot\ and /3 2 {a.i 
and Pi) V-O-V angles disappears , see point (i) above, 
whereas (ii) is still valid, (see Fig. [5Ji); as a result, P 
is suppressed, but is still non-zero. In other words, the 
tT44- spin ordering, imposed onto the centric IAi/amd 
space group, gives rise to an electronic instability which 
ultimately results in a V-V dimerization and formation 
of short and long V-0 bonds, compatible with a stag- 
gered xz,yz orbital ordering. This electronic instability is 
already evident before ionic optimization: the two t {\) 

V sites are inequivalent by symmetry, and, upon ionic 
relaxations, this inequivalency eventually drives the V- 

V dimerization. Also the two oxygens bonded to |,t V 
(or 1,1 V) become inequivalent, in turn giving rise, upon 
ionic relaxations, to a weakly staggered orbital ordering. 
Both effects cooperate to stabilize the polarization. 

Thus the theoretical results confirm, first of all, the 
structure proposed in 12], with dimerization in V-V 
chains in xz— and yz-directions. Furthermore, they show 
that CdV2 04 is ferroelectric. The calculated polarization 
is along the c-direction, and its value is P w 200 /iC/m 2 . 
This is typical for multiferroics with ferroelectricity in- 
duced b y m agnetostriction, like HoMn03 or YMn 2 05, 
see e.g. [14J. The fact that the experimentally observed 
value is smaller, P k, 5 /xC/m 2 , see above, is rather 
common. The same happens in HoMnOa: the measured 
value in polycrystalline material P ~ 90 /iC/m 2 14{ is 
smaller than the theoretical value [IH - which however 
agrees with the estimates obtained from optical studies 
[29J . Most likely it is connected with the poly-crystalline 
nature of the samples: first, polarization is averaged over 



all directions, but most important is that probably in 
these granular materials one does not reach full domain 
orientation during the poling procedure used. 

Why the polarization is absent in the measured sam- 
ples of ZnV 2 04 and MgV2 04, is not completely clear. 
Theoretically we could expect that these spinels, with 
very similar structure, could be also multiferroic. One 
of the reasons could be that these materials have much 
smaller gaps than CdV204 and are close to the localized- 
itinerant crossover 16\. The finite conductivit y o f these 
samples (10 4 times higher than in CdV204 [12j) may 
scramble the results of the measurements. The quality of 
the samples may also matter. In any case, this question 
deserves further study. 

Summarizing, we discovered the first multifer- 
roic ternary spinel with collinear magnetic structure, 
CdV2 04. We thus show that not only a spiral magnetic 
structure can produce ferroelectricity in spinels; the mag- 
netostriction mechanism can do it as well. In the latter 
case, polarization is usually even larger than in spiral 
magnets, which is also the case here. Our study also 
clarifies the very controversial question about orbital or- 
dering and structural distortions in V-spinels, confirm- 
ing that there should appear a strong V-V dimerization, 
that is also responsible for polarization. We conclude by 
suggesting that similar phenomena might occur in other 
spinels, thus broadening the class of multiferroic systems 
to this important group of materials. 

The research leading to the ab-initio results has re- 
ceived funding from the European Research Council 
under the EU 7th Framework Programme (FP7/2007- 
2013)/ERC Grant Agreement No.203523. The exper- 
imental work has been funded by the DFG through 
SFB608 (Cologne). 
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